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ABSTRACT 

Analytical solutions for the three-dimensional Inhomogeneous wave equation 
with flow In a hardwall rectangular wind tunnel and In the free field are 
presented for a stationary monopole noise source. Dipole noise sources are 
calculated by combining two monopoles 180° out of phase. Numerical 
calculations for the modal content, spectral response and directivity for both 
monopole and dipole sources are presented. In addition, the effect of tunnel 
alterations, such as the addition of a mounting plate, on the tunnels 
reverberant response are considered, In the frequency range of practical 
Importance for the turboprop response, Important features of the free field 
directivity can be approximated In a hardwall wind tunnel with flow If the 
major lobe of the noise source Is not directed upstream. However, for an 
omnidirectional source, such as a monopole, the hardwall wind tunnel and free 
field response will not comparable. 

INTRODUCTION 

The relatively high fuel economy available from propeller-driven aircraft 
has renewed Interest In high speed, highly loaded multiple blade turboprop 
propulsion systems. The undesirable features of community noise and, more 
importantly, the high Intensity cabin noise associated with the propellers 
supersonic helical tip speeds have stimulated new theoretical and experimental 
research on the acoustic characteristics of turboprops. 

The acoustic testing of propellers at realistic Inflow Mach numbers can be 
carried out In flight using a suitably scaled model, or performed In a wind 
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tunnel capable of producing high subsonic flow velocities, again using a scale 
model. Since the cost of flight testing Is high, acoustic testing with 
supersonic helical tip speed propellers has Initially been carried out In the 
8 by 6 ft transonic wind tunnel at the Lewis Research Center (1-61 . Figure 1 
displays typical propeller models mounted In the Lewis 8 by 6 wind tunnel. 

The propeller Is driven at the desired speed by an air drive. Noise 
measurements are made with pressure transducers Installed flush with the 
tunnel walls, through the bleed holes shown In the photograph or on added 
structures such as the boundary layer refraction plate. As seen In Fig. 1, 
the hard tunnel wall, wings, or refraction plates could produce significant 
reverberation effects; thus, the accuracy of the noise measurements has been 
questioned. In measurements of this type, the test site Is essentially a duct 
containing a noise source represented by the propeller. It Is not clear that 
this test environment will generally produce results for directivity or 
amplitude which have any relationship to flight test results, although some 
Initial flight to wind tunnel comparisons with the NASA Oryden Jet star 
aircraft ( 7 ,B ) have Indicated reasonable agreement. 

Olttmar (£) has qualitatively argued that flow convection and the highly 
directional nature of a typical propeller noise source will minimize 
reverberation effects when the tunnel Mach number Is greater than 0.6. 

However, Eversman and Baumelster (TO) found the acoustic radiation field to be 
substantially altered by the presence of the duct except perhaps very near the 
propeller. They found a strong resemblance between the acoustic directivity 
on the duct wall and the general shape of the free field directivity for side 
line angles of 15° to 135°. For side line angles less than 45° or greater 
than 135°, similarity between duct directivity and free field could not be • 
expected . 
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Because of computer storage limitations, the results of Ref. 10 were baseo: 
on axial Mach numbers less than -0.5 and source frequencies approximately 25 
percent of the turboprop test frequency (1000 Hz). Considerations In this 
paper will be given to the actual rectangular form of the B by 6 Lewis wind 
tunnel, the design frequency (1000 Hz) and the design Mach number of 0.8. 

It Is the purpose of the present Investigation to mathematically model In 
closed form the radiation patterns of stationary monopoles and dipoles In a 
hard wall version of the NASA 8 by 6 wind tunnel, In free space, and In a half 
space adjacent to a solid Infinite plane. In contrast to the cylindrical 
geometry and low frequency results of Ref. 10, the acoustic response for the 
NASA 8 by 6 wind tunnel will be determined for a rectangular geometry In the 
frequency range associated with the first harmonic of the propeller tests 
(1000 Hz). 

The analysis to follow models the pressure solutions for a stationary 
monopole in an Infinite duct with flow. First, the appropriate governing wave 
equations and boundary conditions will be presented and solved. Next, the 
free field pressure solutions for a monopole source as given by Goldstein (11 ) 
Is developed in terms of the duct coordinate system. These solutions are used 
to describe the resonance frequencies and modal amplitudes of the 8 by 6 wind 
tunnel. Finally, for a monopole and dipole source In the tunnel, the spectral 
response, directionality, and spectral response changes due to. the addition of 
panels or other superstructure Into the tunnel are determined. 

NOMENCLATURE 

A constant of Integration 

a Integer Index, see Eq. (A51) 

B constant of Integration 

b Integer Index, see Eq. (A51) 

C constant of Integration 
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dimensionless sonic velocity, c*/H*f* Eq. (4) 
dipole strength 
frequency, Hz 
tunnel height, m 

VT 

defined Eq. (B) 
wave number, u/c 

wave number In respective coordinates 
dimensionless tunnel width, L*/H* 

Mach number, U*/c* 

Mach number V/c 

source Mach number 

transverse "z" mode number 

number of m modes 

number of n modes 

transverse mode number 

dimensionless acoustic pressure, p*/p*c*2 

time Independent monopole source strength 

dimensionless monopole strength, (H*/c*)q* 

spatially Independent monopole source strength 

equivalent radius, defined Eq. (21) 

radius 

source radius, defined Eq. (24) 
period, sec 

dimensionless time. t*/T* 
axial mean velocity, cM » U*/H*f* 
negative Image of axial velocity, - U 
dimensionless acoustic velocity, see Eq. (2) 
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X function of x 

x dimensionless axial coordinate, x*/H* 

Y function of y 

y dimensionless transverse coordinate, y*/H* 

Z function of z 

z dimensionless transverse coordinate, z /H* 

a nm mode factor, Eqs. (11) and (12) 

(l m defined Eq. (17) 

3 n defined Eq. (16) 

Y Mach number convection factor, Eq. (13) 

6 delta function 

e source-observer angle, defined Eq. (20) 

dimensionless perturbation velocity potential, tj>*/c *H* 
¥ eigen-function 

w dimensionless angular velocity, 2irf*T* 

Su bscripts 

f fixed coordinates 

g Galilean coordinates 

m transverse "z" mode number 

n transverse "y" mode number 

0 space Independent 

s source position 

+ positive x direction 

negative x direction 
Superscrip ts 

* dimensional quantity 

1 coordinate system Eq. (All) 

" coordinate system, Eq. (A25) 
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WIND TUNNEL MATHEMATICAL MODEL 

The first solution to be presented Is for a stationary monopole source In 
a duct with flow as shown In Fig. 2. The equations and boundary conditions 
describing the sound propagation will now be presented along with their 
solutions. 


The linearized gas-dynamic equations for a rnonopole noise source are given 
by Goldstein ((H), Eq. (1.11)) In vector form. For the rectangular duct in 
the fixed coordinate system of Fig. 2, these equations can be combined to 
obtain a dimensionless wave equation of the following form: 


where 


and 





2 2 
a ip 2M a u> 

32* ' c S V X f 


i ,2 

1 a <£ 



q 


3tp 

U a — 

x aXf 


v a 
y ay f 


a tp 



P 


1 §(£_ u 3je_ 

c at f axf 


(l) 


( 2 ) 


(3) 


The usual notation for velocity potential, acoustic velocity, and acoustic 
pressure are used. These and other symbols are defined in the nomenclature. 
Here, the dimensionless speed of sound c Is defined as 


c* 

c a H*f * 


( 4 ) 


The asterisks denote dimensional quantities. 

For a monopole source, wave Eq. (1) can be rewritten as 


2 2 2 2 
, „2t a tp . a cp . a «p 2M atp 

(1 - M) f + — f + —f - r- — £ 


sx 


f 


ay 
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az 


f 


c ax^atf 


2 

c 2 at 2 


= q 0 (t f )a(x f )6(y f - y g ) 4( z f - 


z s ) 


(5) 
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The solution to Eq. (5) can be more conveniently obtained by first rewriting 
Eq. (5) In a Galilean (x ,y ,z ,t ) system moving with the flow, and 

J 3 *1 W 

then transforming Into a Lorentz coordinate system, as outlined by Horse and 
Ingard ((12), pp. 701 and 722). The standard separation of variable solution 
In the Lorentz system and transform back Into the tunnel coordinate system 
shown In Fig. 2 yields 


P(x,y,z,t) « ■— 


N y f5 N z +5 

E E 

n=0 m»0 




^n®m a nm K nm 


Y [1 - MaK sgn( x - xj] 
nm nm c ' 


X exp (X - x s )Y 2 (o. nm K nn , S gn(x - * s ) - h] 


-1a>t 


where 

u* s 2irf * or w = 2ir 

(the frequency f* Is tied to the definition of t a f*t* In the 
dimensionless representation of u) 


( 6 ) 


(7) 


nm 


V , /nirc \ 2 AmrcY 

1 " \*r) - 


* nm (y» z ) = cos(mry) cos pp) 


’? nrn (y s z s ) ° cos < niry <> cos 


a a 1 when 
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N y a Integer,^ 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 
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N z « Integer 

(irj 


(15) 

°n “ 1 

n p 0 R n 

1 

" 2 

n > 0 

(16) 

“m" 1 

m a 0 fl 

m 

1 

° 2 

m > 0 

(17) 


X 

sgn(x - x s ) - jj- 

- x s' 


(18) 


Details of the derivation are given In Appendix A. 

Equation (6) describes the acoustic pressure In the duct for a harmonic 
(e^ ut = e 12irf * t ) source of strength Q q at position x s » y s » z s * 

As seen In Eqs. (6), (11), and (12), the pressure Is composed of a finite 
number of propagating acoustic modes (a nm a l) and an Infinite number of 
nonpropagating modes (evanescent, a nm » 1). At sufficient distance from the 
source, the evanescent modes may generally be Ignored. For this particular 
analysis, approximately five components of evanescent modes are used In the 
numerical calculations. 

Also, as seen In Eqs. (6) and (8) for a given geometry, Mach number, and 
driving frequency f*(c = c*/H*f*)» K nm will go to zero and P will 
approach Infinity (resonance). Figure 3 shows experimental data (13) taken In 
the NASA Lewis acoustic duct laboratory for a compact combination of monopoles 
(small jet). For the resonance frequencies determined from Eq. (8), the 
pressure exhibits the predicted singular behavior; however, real fluid effects 
such as viscosity and nonllnearltles limit the pressure Jump. In addition, as 
seen by the center microphone curve, If the eigen-function (y,z) goes 
to zero, ¥ will dominate over the singularity. This observation Is 
used In the numerical calculations near a singularity. Also, no calculations 
are Performed directly on a singularity. 
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Finally, for a low frequency where only plane waves propagate, the n » 0 


and m <» 0 pressure mode from Eq. (6) Is for x $ equal to zero 


P(x.y.z.t) . !ja n S9"(»n exp J iim a n(j) ,l e -i M t 
dL (1 - H^) [c(l - H^) J 


( 19 ) 


which Is the classic plane wave propagation equation. The convective wave 
amplification Jump at x equals zero comes directly from the H sgn (x) term. 
Free Field Mathematical Model 

The free field pressure solution for a moving <M S ) monopole source Is 
given by Goldstein (H) In Eq. (1.111) In terms of an effective radius Re 
defined by Eq. (1.106), In terms of the x,y,z,t coordinate system of 
Fig. 2, the free field pressure can be expressed for a fixed source In a 


moving fluid (M ■= - M s ) 


P(x,y,z,t) 


Re 


Qq e~^ ut j -In Mfcos 6 f Ml I 

4ifRe( 1 ♦ H cos e) 2 L c " Re ^ * M cos e >j 

- H(x - x s ) » y<x - x s ) 2 + (1 - M 2 )r 2 


1uRe/ c 

o 


( 20 ) 


( 21 ) 


cos 6 + M => 


- * 5 > 
Re 


1 + M cos 0 b 


4 


(x - x s ) 2 ♦ (1 - M 2 )r 
Re 


2 

s 


rf = (y - y s ) 2 + (z - z s ) 2 

Equations (21) and (24) can be determined by reducing the general vector 
equation In Goldstein; however, these equations are explicitly developed 
Morse and Ingard ((12), pp. 718-719). 


( 22 ) 

(23) 

(24) 


by 
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Plane Wall Mathematical Model 


The solution for a monopolo adjacent to an Infinite flat piano can be 
obtained directly from Eq. (20) by combining with a mirrored Image source 
((H), p. 62) of equal magnitude at location - y , 

NUMERICAL RESULTS 

In the calculations to follow, the analytical equations of the past 
sections will be used to determine tunnel resonance, tunnel frequency 
response, source directionality, modal energy content, tunnel alterations, and 
dipole response. The flow Mach number Is taken to be 0.8 in all calculations 
to be presented. In the pressure response plots to follow, the absolute value 
of the pressure squared Is plotted In log coordinates In the form 

10 log (P 2 1) (25) 

The source strength Q q In Eqs. (6) and (19) was arbitrarily chosen as 8.1 
such that the 10 log P 2 had a value In the range of 30 to 70, The 1 In 
Eq. (25) was added to prevent negative value of the log at low frequencies. 

For all practical frequencies 

10 log (P 2 + 1) « 10 log P 2 (26) 

Duct Resonance Freq uencies 

The resonance frequencies of the 8 by 6 wind tunnel were calculated from 
Eqs. (4) and (8) for various n and m modal numbers when K nm Is set to 
zero. The results of these calculations are displayed In Fig. 4, where a 
solid black line represents a duct resonance frequency. Below a frequency of 
100 Hz discrete resonance values can be seen. The first resonance frequency 
occurs at 40 Hz. However, above 100 Hz the resonance frequencies become very • 
dense and the Individual modes are not distinguishable. 

Figure 5 zooms In on the narrow band of frequencies about 1000 Hz where 
many experiments were conducted (1-5) . Since the modes are dense, It would be 
difficult to not excite many adjacent modes with turboprop noise. The BPF 
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(Blade Passing Frequency) tone would not be narrow enough to avoid the duct 
resonances of Fig. 5. 

S pectral Response 

Ihe spectral response of a monopole located In the center of the B by 6 
wind tunnel (x » 0, y o 0.5, z *» 0.375) Is shown In Fig. 6. The 
response Is measured at a position In the corner of the tunnel (x » 0.5, 
y n 0, z » 0), where all the duct modes are present and cut-off modes have a 
chance to decay. Recall the discussion of Fig. 3, by placing a microphone In 
a corner all the resonance modes will be observed. In the frequency range of 
Fig. 6, 1114 duct resonance values occur. 

In Fig. 6, the response of the same strength monopole In the free field 
and adjacent to an Infinite flat plane are also shown. Clearly, the 
theoretical response In the duct of the monopoly Is considerably higher than 
In the free field. Figure 7 displays the same Information In a narrow band 
about 1000 Hz. At low frequencies and certain key frequencies as seen In Fig. 
6, the pressure response In the duct drops below the free field response. 

The V m <y s ,z s ) term In Eqs. (6) and (10) partly accounts for this 
drop. For the present example (Figs. 6 and 7), y g Is 0.5 and z $ /L Is 
0.375; consequently, 

WW - cos (f) cos (f) (21) 

For all propagating modes when n or m Is odd, the value of ¥ will 
be Identical to zero. As a result, the spectral response strongly depends on 
the position of the source. 

So urce Directionality 

Figures 8 to 10 show directionality results for source frequencies of 40, 
90, and 1000 Hz. In Fig. 8, at the frequency of 40 Hz only the plane wave 
propagates. Ahead of the source (x^ < 0) the flow raises the pressure level 
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to nearly 30 while the pressure level trailing the source falls to about 10. 
Thus, a convective amplification of about 20 dB Is seen. The pressure 
variation In the vicinity of the source Is due solely to cut-off modes. 

For the calculation displaced In Fig. 9, four higher order propagating 
modes exist at a slightly higher driving frequency of 90 Hz. The pressure 
near x » 0 Is only a rough approximation In this range. Figure 10 displays 
similar results for a multimodal 1000 Hz source. 

Similar to the propeller In a wind tunnel reported by Fversman and 
Baumelster (10) discussed earlier, for side line angles less than 45° or 
greater than 135° the similarity between duct directivity and free field can 
not be expected. In the free field, the pressure falls off as one over r 
while In the duct the propagating mode remains at constant amplitude. 

Modal Amplitudes 

For a 1000 Hz source, Fig. 11 displays the modal amplitudes as a function 
of their cut-off fi*£Q.’*ncy. As discussed in the Introduction, this closed 
form analyses was to be used to verify three-dimensional numerical approaches 
to the noise source problem which could handle the added complications of soft 
walls and variable area sheared flows. However, the use of conventional 
finite difference or finite element theories would be difficult because the 
modal amplitudes Increase with frequency. Resolution of the near cut-off 
modes would require a prohibitively large number of transverse grid points or 
elements to resolve the harmonic pressure oscillations associated with the 
high cut off frequencies. 

Tunnel Alterations 

As shown In Fig. 1, the wind tunnel is often modified by the addition of 
wall plates or wings, which can change the reverberation characteristics of a 
tunnel. Figure 12 shows the effect of lowering the tunnel wall 0.1 and 0.2 
dimensional increments. As seen In Fig. 12, the average directionality Is 
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flmllar although the level shifts on average by as much as 10 downstream of 
the source. 

Dipole Response 

A propeller, as consldereo In Ref. 10, has Its major noise lobe pointed 
downstream at an oblique angle and a small residual component pointed 
upstream. As suggested by Dlttmar (9). the strong axial directivity of the 
propeller noise source could reduce reverbcraMon effects In the vicinity of 
the source. The monopole just discussed, however, Is an omnidirectional 
source. To eliminate the upstream component of the monopole, a dipole was 
chosen with the alignment shown In Fig. 13, In contrast to the monopole, a 
dipole noise source has strong directionality depending on Its orientation. 

In the case shown In Fig. 13, the pressure lobes are directed toward the side 
walls, and In principle, could be swept downstream. 

The pressure response of a dipole In a tunnel was modeled as two point 
harmonic monopoles 100° out uf phase separated by a small distance 2h. The 
strength of a free field monopole Q q to simulate a dipole of strength D, 
can be expressed as ((H). p. 572) Eq. (13.482)) 



Figure 13 displays the frequency response of this dipole In the narrow 
band about 1000 Hz. As seen In this spectral plot, many of the harmonic peaks 
of Fig. 7 have been cancelled. Figure 14 shows the directionality of a source 
at 1000 Hz, In the vicinity of the source, the dipoles average pressure 
response roughly follows the response of the dipole In the free field. Again, 
as with the propeller analysis of Ref. 10, for side line angles less than 45° 
or greater than 135° similarity between duct directivity and free field could 
not be expected. The amplitudes of the tunnel and free field response are now 
within 10. In an actual transonic tunnel, as shown In Fig. 1, the walls have 
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holes and will therefore appear soft In an acoustic sense. Consequently, tne 
difference between the free field and tunnel results shown In Fig. 14 would be 
the maximum possible. In the actual tunnel, soft wall absorption will occur 
and the tunnel and free field values will approach each other more closely. 
CONCLUDING REMARKS 

Analytical solutions for the three-dimensional Inhomogeneous wave equation 
with flow In a hard wall wind tunnel and In the free field are presented for a 
monopole noise source. Numerical calculation for the modal content, spectral 
response, and directivity for both monopole and dipole sources are presented. 
In the frequency range of practical Importance to turboprop response, some 
Important features of the free field directivity can be approximated In a hard 
wall wind tunnel with flow If the major lobe of the noise source Is not 
directed upstream. However, for a omnidirectional source, such as a monopole, 
the wind tunnel and free field response will not be comparable. 
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APPENDIX A 

Mind Tunnel Solution 

As shown In the body of the report, for a monopole source, the wave 
Eq. (1) can be written as 


n m 2» a£t£ a£<£ 2M _aj > L. 

0 " M> a * 2 ay 2 at 2 ' « a V‘ ' c 2 at 2 

■ q 0 ( t f )6(x f ) A<y f - y s ) A( z f - z s ) 

The solution to Eq. (Al) can be more conveniently obtained by rewriting 

Eq. (Al) In a Galilean (x ,y ,z ,t ) system moving with the flow. 

y y y « 

Therefore, assume 


(Al) 


x, = x + Ut 

f g g 


y f - V 


x + CM t 

g g 


z f - V 


x = x £ - Ut c 

g f r 

t f n t 

f g 


(A2) 
( A3-5) 


As discussed by Horse and Ingard ((H), p. 701), the space derivatives are 
Invariant under this transformation; however, the time derivatives become 




at, 


at ax 

a lie a + fl H 3i£_ „ a«£_ 

at at f ax_ at p at_ ax„ 

g f g f g g 


(A6) 


2 2 

a_s> i_<£ 

atf 


2 2 
o 1 1 - — *P — , 1 12 a JP 

at; a V*g ax 2 


(A7) 


where U equals CM as defined In the list of symbols. Thus, the wave 
Eq. (Al) becomes 

S + S ' ‘'■’‘V 5 '’' * U V* (, « ' V s )S(Z 9 ' V <A8) 

g y g g g 

Finally, to match the general form of equations as used by Goldstein ( ( TjO) , 
p. 45) and Horse and Ingard ((H). P- 721) which employ moving source, let 
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U - ~ V { A9) 

Therefore Eq. (12) becomes 

*4 ♦ *4 * *4 - h *4 > Ut lill - Vt )4<y - y )i(l - l ) (A10) 

ax ay £ az c* ar 0 9 9 9 s g s 

g y g g g 

To obtain a solution of Eq. (A10), Horse and Ingard ((12), p. 722) suggest 
employing a Lorentz transformation of the form 


X' . Y(x g - Vt g ) 


z ' « z. 


Y V 9 “ c 2 \ 


(All-13) 

(A14) 


or combining Eqs. (All) and (A14) 


- Y f ♦ 


(A15) 


where 


■^1 - H 2 


( AT 6 ) 


(A17) 

In contrast to the Galilean transform, both space and time are Involved In the 


transform. Noting that <p[x (.:' ,t') ,t (x' ,t')3 
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a<p( x 1 . t 1 ) a<p ax 1 a<p 

- . . 0 n . . I «... T - 1. I 
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ax' ax„ T at' ax. 


3tp V 3 ip 
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7 - Y 


2 2 2 2 2 

3_j£ _ 2 a. fp 2 a. < p , 2 L 
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^{-X.'-x.t 1 ) 3x1 v 3ifi_ yii 

at " at' at_ ax' at_ “ Y at ' " Y au ' 


at' at ax' at 
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a ip 
ax' 


iLffi = 3 J>_ _ 2y 2 v 


at! 


at 1 


at 1 ax' 


+ Y v 


2,(2 a 2 * 


ax 


,2 


(A18) 
(A19) 
( A20) 
( A21 ) 


16 


«*g - vy - 6 


(f) 


and using Eg. (A15) 


.“a 1 " q o [»' * jj *]) 


( A22 ) 


(A23) 


Thus, the wave equation In the Lorentz system becomes 


£ 

ax 


,2 


ay 1 


az 


2 , Ji 

£_ 

' 2 c 6 at ' ‘ 


<L_!£_ + + Ljsl. _ I_ 9_je_ 3 yq ( yt ')6(x*)A(y' 

* -44. I £ O 


yj)«(z' 


*s> 


(A24) 


2 

The V/c x' term In q Q was dropped because the right hand side of 
Eq. (A24) has a nonzero value only when x 1 equals zero because of 4(x'). 

Finally, the wave equation reduces to an equivalent stationary system when 
It Is assumed 


As a result, 


also 


x' = x"/y 

X" » yx 1 

(A25) 

y" = y‘/y 

y" =* yy' 

(A26) 

z' » z"/y 

z" 3 yZ' 

(A27) 

t* » t"/y 

t" » yt’ 

( A28) 

y 2 (same for y, z, and t) 

ax"* 

( A29) 



( A30) 


(A31) 


(A32) 


17 


Thus, the wave equation In stationary coordinates becomes 


.e - ,c 

3_g> . 3_jp.„ a m i_ a 4> 

A “ A ' A ^ A At 


ax" 2 ay" 2 az" 


?- 22 " rq 0 <t M )4(x")«(y" - y")a(z" - z"> 
z c z 3t „z 0 s s 


Lets now assume the driving source to be harmonic and of the form 


( A33) 


<p(x\y",t") > 4(x",y\z")e 


-lut' 


and 


q 0 (t) » Q 0 e 


-lot' 


( A34) 


Consequently, the governing wave equation becomes 
2 2 2 

+ = y 2 Q 0 *(x")i(y»)4(z") (A35) 

ax"^ ay"^ az"* 0 

where k equals u/c 

The solution to a form of Eq. (A35) Is considered by Horse and Ingard 
((12), p. 500). Because of certain differences, some details of the solution 
to Eq. (A35) are now presented. Assuming 

i * X(x")Y(y")Z(z") 

and substitution Into homogeneous form of Eq. (A35) 


1 alL. . I a ?L_ . k 2 1 a!x_ ,2 

Y ay" 2 2 az" 2 " " * ax" 2 " * 


where k^ Is the separation constant. The solution for X(x") yield 

-Ik x" Ik x" 

X = c^ e x +c 2 e x 


( A36) 


For the Imposed exit condition of no reflection at + and - Infinity, 

+1k x" 

X = c. e x x" > o 

T 

-Ik x" 

• X « c e * x" < o 


( A37 ) 
( A38) 
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Similarly 


which yields 


and 


iai tk 2 k 2 B 1 A_ k 2 

2 " K xY 9 "u 

4 3z „4 X Y a y l( <: y 


Y a A. cos k. y" + A„ sin k y" 
i y d y 


Z a cos k z z" + B,, sin k z z u z" 


( A39) 


<A40) 


(A41) 


where 


,2 ,2 .2 .2 
k z a k “ k x “ k y 


(A42) 

The transverse boundary condition at the hard wall requires the acoustic 
velocities to be zero. According to Eq. (2) at y^ = y" = o and z fc = z" = 0, 
3y/ay " = 0 and 3z/3z" » 0, thus 

$ = C. cos k„ 

+ y 


+lk x" 



cos k z z"e 

x“ > 0 

( A43) 

flk x" 



cos k z z"e 

x" < 0 

(A44) 


The second boundary condition requires the velocity to be zero at y^. «= l 
and = L as shown In Fig. 2. In the coordinate transform, this requires 
that the velocity be zero at y" = y and z" = yL. Consequently, the 


eigen-values k and k are such that Eqs. (A43) and (A44) become 

J *■ 


n m +1k X" 

* nfn (x 11 ,y",z") = C + cos r COS z" e A 


-Ik x" 

» C cos — y" cos ~ z“ e x 
Y YL 


X" > 0 


X" < 0 


(A45) 
( A46) 


where 


- (?) - (?f ) 2 - 


Y«o J J 


(A47) 
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Equations (A45) and (A46) represent the separate homogeneous solution to 
the wave Eq. (A35) In the negative and positive domain of the spatial 
variable x. What remains Is to determine the constants C f and C_ from 
the nonhomogenous term on the right hand side of Eq, (A35). 

For the no flow case considered here, the x momentum equation can be 
Integrated between x ■* - c and x ** + c. As e goes to zero, It Is seen 
that both the pressure and the velocity potential are continuous across a 
monopole. Consequently and C are equal in magnitude and 
Eqs. (A45) and (A46) become 


$ 

nm 


(x",y" 


z") 


nm 


cos 




Ik lx" | 
e 


(A48) 


In contrast, a dipole would have continuous velocity and a discontinuous 


pressure and velocity potential. 

The last step remaining Is to determine the constant C nm . For 
simplicity, let 

1k y |x"l 

* o c e * 

Y nm nm 


( A49) 


Thus, 


00 C3 

1 - Z c ° s 


nir 


niff 


y" cos r z' 
Y yL 


( A50) 


n=o m-o 

Substituting Eq. (A50) Into the wave Eq. (A33), multiplying both sides of the 
equation by cos airy/y cos birz/yL and Integrating yields for any value of 
norm 
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The above Integration will remove the delta function from the differential 
equation, because for any function g 

f 4(y" - yj)g(y")dy - g(yj) 

Finally, because of the orthogonality of the cosine, Eq. (A51) becomes 


( A52 ) 


* [* 2 - ( f ) 2 - ( S ) 2 ] ♦» 


Q 0 A ( x " ) cos y» cos z" 


where 


n r, - 1 


n <* o 


On ■ t 


fj ■ 1 m « 0 
m 


(3, 


ra 2 


Y J $ 

L GA 

n m 


n > o 


m > o 


y I. ( A 53 ) 


(A54) 

(A55) 


The multiple values of 13 n are a consequence of the properties of the 
2 

Integral of cos ny". 

Finally, substituting Eq. (A49) back Into Eq. (A53), Integrating both 
sides of the equation between x «* - e and x « + c and letting c approach 
zero (see (12). p. 133 for details) yields 


'nm 


1 Q n cos f y» cos jjj* Z ; 

2 Wm 


(A56) 


For ease of calculation let k x In Eq. (A47) be rewritten as 


k o — K a 
x c nm nm 


( A5 7 ) 


as defined by Eqs. (8), (11), and (12) In the body of the report. Thus, 

oa_ co 

4(x",y",z") - 22 22 

ll 


n=o m=o 


IQ, 


J nir „ 

cos r y « 


2wa K L(3 |3 
nm nm n m 


mtr . 
cos t z' 

xL-.ll 


ntfv 11 nurz 11 D la nm K nm c ^ 

Y cos yL 6 


( A 58 ) 
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Equation ( A58) represents the potential In the stationary system. This 
equation must be transformed back Into the tunnel coordinate system In 
Fig. 2. First, the equation Is transformed back Into the Lorentz's system by 
substituting Eqs. (A25) to (A28) Into Eq. (A58) to yield 


CO qo 


mirz^l 


-IQ c \ 1 cos mry£ cos — — cos mry' cos *-j* 


rntrz 1 


n m nm nm 


nao m»o 


X 1~ Y I X ' |tt m k nrn /AKQ\ 

„ c ■ 1 ’ nm nm ( A59) 


where Eq. (A34) 


q>(x' ,y' , z 1 ,t') « ^(x'.y'.z'Je -1 ^ 


Next, Eqs. (A59) and (A60) must be transformed to the Galilean 
( x ,y , z ,t ) system. Multiplying Eq. < AT 4 ) by y 

y y s u 

2 V 2 

■ Ai - ? V 


(A60) 


(A61) 


and noting that 

2 2 l q 2 

t 3 . t s (l - H 2 ) . H 2 t g = -f . H t, = ? t g 

then replacing t In Eq. (A61) with Eq. (A62) 
y 

T*' ■ *g - < X 8 - V V 

Substituting Eqs. (All), (AT 2) , (A3), (A63), and (A9) Into Eqs. (A59) and 
( A60) yields 


( A62) 


( A63) 


DO GO 


‘<YV° ■ 


-ly “rnnW l* 9 * V 0 t g | 

2L “ WinJ ^nVnm a nm 
n=o m=o 


( A64 ) 
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-i„t„ uvW <x g * U 0 t g ) 

*<VWV " ®^ x g ,y g' z g^ 0 9 e 


( A65) 


where 


rmrz 

W z s> " cos n «»s cos — 


{ A66) 


mirz. 


*nm (y g’V * COS niry cos -j-a 


( A67 ) 


The sign of (x + U Q t) In Eq. (A65) need not be considered because this will 
be removed In the next transform to tunnel coordinates. 

Substituting Eqs. (A2) to (A5) Into Eqs. ( A64) and (A65) yields 

«„ y , , . -l5fic „ ’c “nn. K nn,' x fl (A68) 

* (x f' y f z f> 2 U LtLj 


n>=o m=o 


n’ m nrnm 


~1u>t ( 


1o>y 2 V 


*p( ^p»y^r Zjrtt|;) o 4(Xj;,y|;,Zj:)e B 


2 A f 


( A69) 


The w terms have the same functional form. Combining Eqs. (A68) and (A69) 
and substituting Into Eq, (3) In the body of the report yield for pressure 


P(^j;>yj:»2j:tt^) 


Ul UJ 

5a T VjteM 2 r, 

2L LjLj T V ' 

‘"[/V ‘“JWV ' Hx f>] 

IxJ 


H a„ m K„ m sgn (x^) 1 
nm nm 3 f J 


n=o m=o 


Note, In taking the spatial derivative of e 


. lut. 


A|X f' 

ax e 


= sgn (x f )A e 


A I x | 


where 


(A70) 


(A71) 


sgr,(x f ) =7^ 


<A72) 
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as easily seen by rewriting two sepatate equations for exp A|x| for x 
greater than zero and x less than zero. Finally, we transform to the tunnel 
coordinate systenw noting that 

x « x f + x s Xj. « x - x $ (A73) 

y «> y f , z - z f , to t f ( A74-A76) 

The resulting equation Is Eq. (6) In the body of the report. 


24 



REFERENCES 

1. Dlttmar, J.H., Blaha, B.J., and Jerackl, R.J., "Tone Noise of Three 
Supersonic Helical Tip Speed Propellers In a Wind Tunnel at 0.8 Mach 
Number," NASA TM-79046 (197B). 

2. Dlttmar, J.H., Jerackl, R.J., and Blaha, B.J., "Tone Noise of Three 
Supersonic Helical Tip Speed Propellers In a Wind Tunnel," NASA TM-79167 
(1979). 

3. Dlttmar, J.H., Jerackl, R.J., "Additional Noise Data on the SR-3 
Propeller," NASA TM-81736 (1901). 

4. Dlttmar, J.H., Stefko, G.L., and Jerackl, R.J., "Noise of the 10-Bladed, 
40° Swept SR-6 Propeller In a Wind Tunnel," NASA TM-82950 (1982). 

5. Dlttmar, J.H., Stefko, G.L., and Jerackl, R.J., "Noise of the 10-Bladed, 
60° Swept SR- 5 Propeller In a Wind Tunnel," NASA TM-83054 (1983). 

6. Dlttmar, J.H. and Jerackl, R.J., "Noise of the SR-3 Propeller Model at 2 
and 4° Angle of Attack," NASA TM-82738 (1981). 

7. Mackall, K.G., Lasagna. P.L., Walsh, K., and Dlttmar, J.H., "In-Flight 
Acoustic Results from an Advanced-Design Propeller at Mach Numbers to 
0.8," AIAA Paper No. 82-1120 (1982). 

8. Dlttmar, J.H. and Lasagna, P.L., "A Preliminary Comparison Between the 
SR-3 Propeller Noise In Flight and In a Wind Tunnel," NASA TM-B2805 
(1982). 

9. Dlttmar, J.H., "Why Credible Propeller Noise Measurements are Possible In 
the Acoustically Untreated NASA Lewis 8 by 6 Ft Wind Tunnel," J. Acoust. 
Soc. Am., vol. 75. no. 6, June 1984, pp. 1913-1914. 

10. Eversman, W. and Baumelster, K.J., "Modeling of Wind Tunnel Wall Effects 
on the Radiation Characteristics of Acoustic Sources," AIAA Paper 

No. 84-2364 (1984). 


25 



11. Goldstein, H.E., "Aeroacoustlcs," McGraw-Hill, New York, 1976. 

12. Horse, P.H. and Ingard, K.U., "Theoretical Acoustics," McGraw-Hill, New 
York, 1968. 

13. Sued, G.P., Gaumelster, K.J., and Ingard, K.U., "Interaction of a 
Turbulent-Jet Noise Source with Transverse Modes In a Rectangular Duct," 
NASA TP-1248, June 1978. 

14. Reynolds, Douglas D., "Engineering Principles of Acoustics Noise and 
Vibration Control," Allyn and Bacon, Inc. (1981). 


26 



OIWttklAL PACK. .S 
OF POOR QUALITY 



PROPELLER ALONE AND AT ANGLE 
ATTACK 



WING SHIELDING 



INSTALLATION EFFECTS AND BLADE 
MOUNTED TRANSDUCERS 
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Figure 1. - Advanced prop-fan installations in NASA Lewis Research 
Center 8 x 6ft Transonic Wind Tunnel. 
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Figure 3. - Modal response ot a small rectangular duct (rel. 13). 
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Figure 4. - Resource frequencies of the 8x6 ft wind tunnel. 
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Figure 5. - Resonance frequencies of the 8x6 ft wind 
tunnel in a narrow band about 1000 cycles per sec. 
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Figure 6. - Spectral response of a monopole located in the 
center of the tunnel (x$ • 0, y s * 0. 5. z s * 0. 375) at 
x ■ 0. 5, y ■ 0. z * 0. tO • 8. 1, M ■ 0. 8) 
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Figure 7. - Spectral response of a monopole of figure 6 
in a narrow band about 1000 Hz. 
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Figure 10. - Wall pressure response for a 1000 Hz mono- 
pole (y “ 0, z » 0.375, x s B 0, y s - 0.5 and z s * 0.375). 



Figure 11. - Modal amplitude of propagating acoustic modes 
for 1000 Hz monopole source. 
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Figure 12. - Effect of tunnel alterations on directivity of 
1000 Hz monopole source (x s ° 0.0, y s = 0.5, z s = 0.375, 
y - 0. 0, z - 0.375). 
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Figure 13. - Spectral response of a dipole in a narrow band 
about 1000 Hz (Q = 8.1, 2h - 0.05, x B 0.5, y - 0, z - 0, 

M » 0. 8) 
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Figure 14. - Wall pressure response for a 1000 Hz dipole 
at y - 0 and z ° 0 (x s B 0, y s c 0. 5, z s = 0. 375, M - 0. 8), 
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